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An analytical study is presented of the crossover in the gapless attractive XXZ 
chain from antiferromagnetic to ferromagnetic behaviour at low to high temper- 
ature, respectively. In particular, an analytic formula for the crossover in the 
long distance asymptotics and explicit results for the nearest-neighbour longitu- 
dinal correlation are obtained. We also provide results for the specific heat and 
magnetic susceptibility for various anisotropies. 



1 Introduction 

In this paper we pursue the question to which extent the partially anisotropic 
Heisenberg chain with nearest-neighbour coupling 

1 L 

h xxz = 5 E ( CT f CT f+i + "Y°T+i + A(j f (!) 



known that ([j]) for anisotropy 
however with quite different 



is described by conformal field theoryEJti It 
— 1 < A < +1 has a critical groundstateB 
physics regarding the elementary excitations E3tij. There are gapless elemen 
tary spin excitations throughout — 1 < A < +1, but additional bound states 
in the "attractive" parameter regime — 1 < A < 0. Note that A = — 1 cor- 
responds to the isotropic ferromagnetic Heisenberg chain with fully polarized 
groundstate. 

Correlation functions are difficult to calculate, ^Lnalytical results are known 
for the nearest-neighbour longitudinal correlations (cr? <t? +1 )t=o (which will 
be generalized to finite T in Sec. |]) as well as fkjld-theoretical results for 
long-range behaviour of correlations for T = OBQ and small non-zero T 
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Amazingly, the nearest-neighbour correlation (<jf<jf + i)T=o is always negative 
for all anisotropics — 1 < A < +1. In addition, the longitudinal correlation 
(cTq a^)T=o for large distances r is negative in the attractive regime. 

The finite -temperature results for (a^a^) as found in numerical studies 
are quite rich 113. For fixed separation r / a transition^ from negative to 
positive values appears for some temperature Xb(A,r). InE3 the question was 
raised whether To (A, r) possesses a well defined (and non-zero) value in the 
limit t — > oo. Alternatively, in an asymptotic expansion in exponentials 

(a V)=^ex P (-r/0 + ..., (2) 

the temperature dependence of the amplitude A of the dominant term was 
addressed and it was argued that there is a sign change for some well defined 
temperature Tb(A) (— Tq(A, oo)). Nothing like this is observed in the "repul- 
sive" regime where correlations show simple antiferromagnetic oscillations. 

In the current work we want to understand analytically the sign change 
phenomenon. In particular, we aim at an analytic formula for To, see (p8|). 
For two reasons, these questions are not purely academic. First, the "attrac- 
tive" XX Z chain is often considered as a system with ferromagnetic interac- 
tions, however with antiferromagnetic groundstate, i.e. vanishing magnetiza- 
tion. Our analysis is aimed at a resolution of the somewhat unintuitive picture 
in the way of a crossover from dominant ferromagnetic correlations at high to 
dominant antiferromagnetic correlations at low temperatures. Second, we ad- 
dress the fundamental problem of the additional energy scale that has to enter 
the description of the "attractive" regime, but is absent in the "repulsive" case. 
Naturally, we will be led to pay attention to the existence of bound states. 

In Sec. we review the known groundstate properties of the XX Z chain, 
and in Sec. |2] we set up our formalism for finite temperatures on the basis of 
the quantum transfer matrix. In Sec. ^ we address the problem of crossover 
phenomena in correlation functions and macroscopic properties. Finally, in 
Sec. | we summarize our results and list the open problems. 

2 Groundstate Properties 

In (Q) the anisotropy A is conveniently parametrized by 

A = cos 7 (3) 

where the range < 7 < n/2 corresponds to repulsive interactions and 
7r/2 < 7 < 7r to attractive interactions. The elementary excitations on the 
(antiferromagnetic) groundstate are free states with dispersion relation 

e/(fc) = wsin/c, < k < it, (4) 
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and "sound" velocity 

V = TT. (5) 

7 

In addition to these "free states" there occur "bound states" in the attractive 
regime tt/2 < 7 < it. Depending on the anisotropy there are different bound 



states labeled by an integer \x in the range 1 < /j, < 



with dispersion 



e b (k) = u2sin-yi + a2 s in 2 -, a M = cot ( pi— — -J . (6) 

and < k < it. Note the same velocity of the bound state dispersion as for the 
free state dispersion. However, for finite momentum the energy of the bound 
states is larger than that of the free states. 

For the field theoretical description of the XXZ chain the occurence of 
bound states in the attractive regime does not have any fundamental conse- 
quences, because of identical velocities of free and bound states in the long- wave 
limit. Within the bosonization approach the XXZ chain is mapped to a Sine- 
Gordon model where the interaction term is argued to be infrared irrelevant 
(resulting in a Gaussian model). In this way correlation functions for T = 
and small T > with asymptotic behaviour 

are derived Bi. To our knowledge the breakdown of this picture at higher 
temperatures in general, and the occurence of a new temperature scale in the 
attractive regime (with bound states) in particular has not been studied for 
the Sine-Gordon model. 

The scaling dimension x and "lattice momentum" Pq are generally different 
for different correlation functions and given by a formula obtained from lattice 
calculations 

x = l_yA 5 2 + 1 m 2 + k _ P = (S- m)7T, (8) 
2 2(1-7/71-) 

where S, to, k are integers corresponding to spin, lattice momentum and posi- 
tion in the conformal tower. For the longitudinal spin-spin correlation function 
the selection rules enforce S = 0, but leave open the values for to and k. A sim- 
ple inspection shows that (to, k) = (1,0) (= (0,1)) gives the smallest scaling 
dimension x for the repulsive (attractive) regime. 
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3 Finite Temperatures 

3.1 Quantum Transfer Matrix 

For a trpatmpnt nf finitp temperatures we employ a convenient transfer matrix 
approach EJo'ESEj. To this end the quantum chain at finite temperature is 
mapped via a Trotter-Suzuki decomposition onto a two-dimensional classical 
model on a square lattice of width L (= chain length) and height N (= Trotter 
number) and staggered interactions. The free energy and the decay of static 
correlation functions is completely described by the largest eigenvalue (InAo) 
and the next-largest eigenvalues (lnA^) of the column-to-column transfer ma- 
trix 

/ = --jlimlnAo 

c - = *(&) r+ *(£)' + ™ (9) 

In the limit N — > oo this matrix is referred to as the quantum transfer matrix 
(QTM) being the closest analogue to the usual transfer matrix of classical spin 
chains. 

For the XX Z chain integrability is manifest at the level of the QTM in the 
following way. There is a commuting family of matrices T{x) generated by the 
spectral parameter x. The QTM is identical to T(0). As most of the physical 
properties are directly given by the logarithm of T(0) we define H = — lnT(O). 
The lowest eigenvalue of 7i gives j3f where / is the free energy per chain 
site. Furthermore, the QTM enjoys translational invariance along the vertical 
axis within the two-dimensional geometry. The corresponding momentum V 
operator is generated by T{x) through differentiation. In summary we have 

H = -lnT(O), 

V = -isin~f-^hiT(x)\ x =o- (10) 

Note the similarity to the integrability structure of the Hamiltonian H and mo- 
mentum operator P with one important difference. The momentum operator 
P is given by the row-to-row transfer matrix and H is given by its deriva- 
tive. Consequentially, we will often observe a formal correspondence TL 1 V 
< — ► iP, iH. Of course, in a suitable continuum limit this is easily understood 
as Hamiltonian and momentum operator are evolution operators with respect 
to orthogonal directions. 

Next, we review the already known eigenvalue equations for the operator 
T(x). The corresponding eigenvalues are denoted by A(x) from which the 
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eigenvalues £ and K, for H and V are directly obtained 
£ = -lnA(Q), 

JC = -isin7-^-lnA(x)| x =o, (11) 
ax 

the latter one taking all integer multiples of 27r//3. This discretization of K, 
naturally entails a discrete spectrum for £ . 

3.2 Eigenvalues and N on- Linear Integral Equations 
The largest eigenvalue is given by@ 

bAW ~«,, + ££j^L*. (i2, 

where eo is the groundstate energy and A = 1 + a, A = I +a are solutions to 
the following set of non-linear integral equations (NLIE) 

^a(x) = ^— + - - Bh 

v ; cosh^x 2(tt-7) 



OO 



+ / \k(x — y) In A(y) — k(x — y — + ae) In -A (2/)] dy. (13) 
J —00 

h is the externally applied magnetic field and the integration kernel is given 
by 

1 f°° sinhff -7) kcos(kx) , . 

k(x) = — / ^ ^ —dk. 14 

v ; 2nJ_ 00 2cosh^fcsinh^pA; V ' 



The corresponding equation for a is obtained from Eq. ([L3|) by exchanging 
i — > — i, h — ► —ft and a, ^4 <-> a, A. These equations will be studied numerically 
in the next section yielding results for the specific heat, magnetic susceptibil- 
ities, and nearest-neighbour correlation for various anisotropies A and wide 
temperature ranges. 

The next-largest eigenvalues describing the decay of longitudinal spin-spin 
correlation-fjjnctions, i.e. eigenvalues with spin quantum number S = 0, are 
given byNlia 



lnA(x) = -f3e (x)+\n 

, 1 f°° HAA(y)}dy 



7T 7T 

tanh — (x — 0i) tanh — (x — 82) 

27 27 



27 Joe cosh^(x-y)' 



(15) 
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where A and A are determined from 

ma x = h 7ri H ; -p/i + m — ; r-rr 

{ > cosher 2(tt- 7 / silm ^^-(yo-n/S)) 

-if (a: - (0i + i 7 /2)) - K{x - (0 2 + if/2)) 

/•OO 

+ / [fc(x — y) lnA(y) — fc(x — y — if + ie) hi A(y)] dy. (16) 

J oo 

and K(x) is defined by K(x)' = 2nik(x), or explicitly 

r°° sinh(^ - 7 ) fcsin(fcx) 

= * / V2 = ' i^^dk. (17) 

J_oo 2fc cosh ^fc sinh — j-^fc 

There appear three parameters 6±, 6* 2 , and j/q m the upper equations corre- 
sponding to "hole positions" (both 61 2 on the real axis, or forming a com- 
plex conjugate pair) and one "complex rapidity" (with Im(yo) = 7r /2) of the 
underlying Bethe ansatz pattern for the QTM. There are higher-lying states 
satisfying similar integral equations, however with more parameters 6i and yj. 

The parameters 9±, 62, and yo do not take arbitrary (continuous) values, 
they have to satisfy the coupled equations 

a(0i + i-y/2) = a(6 2 + if/2) = a(y + if/2) = -1, (18) 

leading to a quantization of the eigenvalues. In general the equations ( |l6| , [i"8| ) 
have to be solved numerically in order to deal with the problem of non-linearity. 
For the case of the largest and next-largest eigenvalues an iterative approach 
proved useful showing convergence within a numerical accuracy of 10~ 6 already 
after approx. 10 steps, see ref. Il3 for the non-critical cases with A > 1 and 
A < — 1, and ref. Illl for the critical case A = 1. A comparable numerical 
analysis for —1 < A < 1 is in progress, however not yet completed. For 
this reason we apply an analytical study within a reasonable approximation to 

Oil- 

3.3 Conformal Field Limit for Low Temperatures 



IntB the equations ( |l5| , [16|) were treated in the low-temperature limit. Despite 
the apparent complexity of the non-linear integral equations, quite a universal 
picture evolved. A certain symmetry of the integration kernel allowed for 
analytic manipulations avoiding the necessity of an explicit solution of the non- 
linear equations. Each eigenvalue could be written in terms of dilogarithmic 
integrals which resulted in the explicit formula 

2ir 

In A = -0ea - —T{x - c/12) + o(T 2 ) + iP , (19) 
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where the central charge is c = 1 and the scaling dimension x and lattice 
momentum Pq are given by (Q). From ([|,|l9]) we see 

C r ~ cos(P r)e ~ xl r . (20) 

which coincides with the conformal field theory result (Q) in the low-temperature 
limit. 

On one hand we see that ( p^| , ^6| ) recover CFT, on the other hand, the 
non-linear integral equations are not restricted to low temperatures, but go 
beyond. It is this behaviour that will be studied in the next sections. 



3-4 Solution in Lowest Order 

A principal treatment of the NLIE is based on iterations. For instance in 
the case of ( |l3| ) a reasonable initial choice for the functions a and a is just 
ciq = ao = 0. Inserting this into the right hand side of ( |l3| ) leads to 

ta "W— sin- < 21 > 

7 

This process should be continued ad infinitum, however we content ourselves 



with the approximation a\ for the function a. Inserting a\ and a\ into (12) 
we obtain the first order approximation (10A) in the sense of an iterative 
procedure to the largest eigenvalue of the QTM which already reproduces 
correctly: (i) the low-temperature asymptotics In A = — (3eo + c(irT/6v) with 
central charge c = 1, and (ii) the high-temperature behaviour In A = In 2. 
This success may be understood in the following way. At low temperatures 
the corrections described by the integral terms in (|l^) are small, and at high 
temperatures both of them cancel each other. We conclude that the (10A) is 
a respectable approximation useful for the entire temperature range. 



Free States 

Next, we study the excitations (|l6|Jl8|) in (10A) yielding 



In A(x) = In 



7T 7T 

tanh — (x — 61) tanh — (x — 82) 
27 27 



(22) 



where we have dropped the common offset — j3eo(x) in (fl5|). From this and 
( |i"T| ) we derive the eigenvalues 



K, = Kf(ei)+Kf(e2), K/(e) = iusinhe, 



(23) 
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where the parameters 9i have been parametrized by e^(> 0) such that Ki is a 
unique function of Ci except for the sign in the n / dispersion which is given by 
the sign of the parameter 8i . Note the similarity of the "momentum-energy" 
dispersion K/(e) for the "free states" of the QTM with the dispersion (Q) 

= ( 24 ) 

We know already from general principles that K. is strictly equal to an integer 
multiple of 2tt//3. This does not necessarily imply a similar property for each 
individual re, = Kf(ej), however values close to multiples of 2tt / (3 are generally 
taken. 

We see this most easily for the case K\ — —K2, i- e - $1 — —82- The quan- 
tization condition ( |l8| ) for yo requires yo = 00. If we apply the (10A) treat- 
ment to (fl6|), the quantization condition ( |l8| ) simply reads h\a{0i +vy/2) — 
i/3v/ sinh (ir/^Oi — odd multiple of iri. Using this when inserting ( p2[ ) into 
(|ll|) we directly find ki = — &2 = 27r//3. 

If both ki and K2 take the smallest possible values of same sign and the 
XX Z chain is in the attractive regime 7r/2 < 7 < tt the quantization conditions 
impose complex values for the parameters #1,2- This leads us to the study of 
bound states of the QTM. 



Bound States 

Here we content ourselves to a motivation of why in the attractive regime 
complex solutions to (f8h close to 



n,2 



= e±i 7- 



(7-5), m 



appear. To this end let us assume that 9\ has a positive imaginary part and we 
calculate lna(#i +ij/2) from (fl6|). For low temperatures the dominant term is 
— (3v/ cosh (ir/j)(6i + ij/2) whose real part tends to +00 in the limit /3 — > 00. 



The only term matching this divergence is the yo term in the first line of (16) 
if 81 — yo = 7 — 7r. Likewise we conclude for 62 with negative imaginary part 
the condition 62 — yo + 7 = 0. Using Imyo = tt/2 we write yo = + in/2 and 
arrive at (|25|). 

We like to note that for the case of higher-lying states involving more than 
one complex rapidity yo there are additional bound states. These always in- 
volve two complex conjugate parameters and a "string" of complex rapidi- 
ties yj with minimum number 1 and maximum number [7/(71" — 7)], a situation 
similar to that of the Hamiltonian. However, for the crossover phenomenon in 
the correlation functions the case explicitly studied above is sufficient. 
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The momentum-energy dispersion is obtained from (|22j,|25|) inserted into 
( pi] ) and eliminating 9. More easily we may use the relation to (^) Kb(e) = 
—ieb(ie). With /i = 1 we find 



Kb(e) = ±v 2 sinh ^ J I — a 2 sinh 2 ^, a = cot ( — — ~7M • (26) 



4 Crossover Phenomena 

In this section we are going to investigate various crossover phenomena in 
correlation functions as well as macroscopic properties of the XXZ chain in 
the attractive regime. Of prime interest is the calculation of the crossover 
temperature as function of the anisotropy parameter — 1 < A < 0. 



4-1 Crossover in Correlation Functions 



With the knowledge of Sec. 3.4 we are in the position to examine the crossover 
of the asymptotics of the longitudinal correlation function from ferromagnetic 
behaviour at high temperature to antiferromagnetic behaviour at low temper- 
ature. In particular we want to derive an analytic formula for the crossover 
temperature as function of the anisotropy A. In principle such a crossover may 
occur under two separate circumstances 

• (i) the prefactor Ai of the dominant term in (||), i.e. that with longest 
correlation length £x = — £o)i turns zero for some temperature. 

• (ii) the leading and next-leading terms in (^|) have prefactors A\ and 
A 2 of different sign, and a crossover of the correlation lengths £i and £2 
occurs for some temperature. 

Unfortunately, we cannot calculate the prefactors-as matrix elements are cur- 
rently out of reach of our formalism. In the workE3 scenario (i) was advocated, 
because it naturally leads to a crossover for (cif^n) & t a temperature T)(n) 
which in dependency on the separation n converges to To with exponential rate 
(formula (2.10) of E2l with three fit parameters). For scenario (ii) the depen- 
dence of T>(n) would be algebraic Tq + A/n with two fit parameters To and 
A. In fact, the exponential fit looks more convincing than the algebraic fit. 
However, the latter procedure involves one fit parameter less. 

As the implications of the fit procedures are not conclusive we performed 
a numerical study of the QTM for finite N which already allows for quite 
accurate, non-trivial high temperature investigations for values of A close to 
0. Surprisingly, this exercise having been aimed at discriminating between the 
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two scenarios showed for N — 4 that actually both take place at the same 
temperature! However, for N — 6 only scenario (ii) is realized. We have 
performed complete eigenvalue computations for much larger numbers N up 
to 14 and always find scenario (ii) realized. For these larger numbers of N 
we have not yet computed the matrix elements such that we cannot decisively 
comment on (i) at the moment. 

In the following we investigate scenario (ii) in detail, i.e. the crossover 
in the two largest correlation lengths, or equivalently in the corresponding 
"energies" £ for the lowest state with K. = and the two lowest states with 
K, = ±27r//3 which are free and bound states, respectively. Such a level crossing 
analysis can of course be performed within our approach. 

The condition for the crossover of a free state with momenta K\ = — K2 = 
2-k/Po and a bound state with momentum k = 2n/f3o where /3 corresponds to 
crossover temperature To is 

e = ex + £2 and 2k = Ki + K2 



=*> 2 sinh e l \Jl-a 2 sinh 2 e 4 = sinhe^, (27) 
with solution sinh 2 ej = Inserting this into t/(ei) = 2n / (3$ yields 



\/3 sin 7 



IT — 7 7T 



T ^x-r ta H-T2j- (28) 

This is in excellent agreement with the numerical analysis of0, see Fig. |[ 
An almost perfect agreement with the numerical results has been achieved by 
introducing a multiplicative correction close to 1 and independent of A. Of 
course, some correction term had to be expected as we only calculated within 
the 10A. 

In addition we show "analytical" results for the nearest-neighbour corre- 
lations (<Tj <Tj +1 )t>o derived as derivatives of the free energy with respect to 
A, see Fig. |. 



4-2 Crossover in Macroscopic Quantities 

Finally, we like to present a qualitative argument based on the dispersion rela- 
tions of the elementary excitations of the Hamiltonian (^,||) to show crossover 
between ferromagnetic and antiferromagnetic behaviour at high and low tem- 
peratures, respectively. Notice that the dispersion of free states (|J) is quasis- 
linear, i.e. is approximated well by a linear relation as long as the momentum 
transfer is less than the reciprocal lattice vector tt. For the bound states (ph 
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Figure 2: Temperature dependence of the nearest-neighbour correlation ( IJ f eT f-^i) f° r A = 
— 1, —0.9, ...,0.9, 1.0. Note that for < A < 1 and all temperatures the values are negative. 
For — 1 < A < there is a crossover from positive to negative values when passing from high 
to low temperatures. For A = — 1 and all temperatures the correlation takes positive values 
which in the limit T = should approach the value +1 in a singular manner. 
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Figure 3: Depiction of x{T) ■ T = J2 ( S S f+l) for ( a ) A = _1 > ~ - 9 ' ->° and ( b ) A = 
0, 0.1, 1. Note the existence of maxima of x(T) ■ T at temperatures of the order Tq for the 
attractive case (a). No maxima are observed in the repulsive case (b). 

the linear regime is much smaller. For < k < (jr — 7)^/27 the dispersion be- 
haves like £f,(fc) ~ i?2sin§, i.e. it is quasi-linear. For (tt — 7)77/27 < k < it we 
find the dependence £(,(&) ~ v 2a sin 2 k/2, i.e. quadratic behaviour typical for 
the elementary excitations of an isotropic ferromagnetic system. The crossover 
takes place for asinfc c /2 = 1 with crossover temperature 

T& = e b (k c ) = 2^ 2 v/a = 2 3 / 2 tt^ tan ^1 . (29) 

Note that Tq is larger than To by a factor of 20.52..., though the functional 
dependence on A is identical. We see that both arguments (crossover of free 
and bound states of the QTM, and crossover from linear to quadratic dispersion 
of the bound states of the Hamiltonian) share the same physical origin, but 
apply to different properties. In fact, at temperatures of the order Tq we see 
characteristics in the susceptibility data which do not exist for the repulsive 
regime < A < 1, see Fig. [|. For any — 1 < A we have lim-r^o xCKlr-T = 0, 
whereas for A = — 1 the result is divergent, x(T) — J/(6T 2 ), seeE3'E2l. The 
noncontinuity of the limiting values is the reason for the observed temperature 
maximum of the quantity xCO ' T '. 

For completeness we show results for the specific heat of the XX Z chain 
in the attractive as well as repulsive regime. Note there are no characteristics 
like finite temperature maxima in c(T) in addition to the usual "peak" typical 
for spin systems with a finite number of degrees per lattice site. 
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Figure 4: Depiction of c(T) for (a) A = -1, -0.9, ...,0 and (b) A = 0,0.1, 1. 

5 Conclusion 

We have studied the competition of ferromagnetic and antiferromagnetic order 
in the longitudinal correlation function of the spin-1/2 XX Z chain in the at- 
tractive regime. Within the scenario of level crossing of the quantum transfer 
matrix we found an explanation of the sign change in the dominant exponen- 
tial asymptotics from positive to negative values when passing from high to 
low temperatures. Notably, we found an analytic formula for the crossover 
temperature To (A) (|28| ) as function of — 1 < A < which may be viewed as 
the maximum temperature below which conformal field theory is applicable. 
The agreement of To as determined numerically and analytically is very good. 

The physical origin underlying the crossover phenomena in correlation 
functions (as well as the magnetic susceptibility) may be viewed in the existence 
of bound states. The various manifestations of these peculiar properties of the 
spectrum of the quantum transfer matrix (as well as the Hamiltonian) have 
been discussed. A lesson to be drawn from these findings is that the physical 
properties of systems possessing bound states are much richer than the usual 
low-energy treatment within bosonization and a subsequent approximation by 
Gaussian models. 

We want to point out two important questions which remained unan- 
swered. First, for the explanation of the crossover there are two different 
scenarios conceivable. For the numerical analysis^ a sign change of the coeffi- 
cient A\ of the leading term in (||) was assumedE3, the analytical work in this 
paper was based on a level crossing scenario as also found in explicit numerical 
treatments of QTMs for Trotter numbers 4 < N < 14. The new level crossing 
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scenario would suggest a numerical size analysis of the crossover temperature 
somewhat different from that employed in Ed. In preliminary numerical work 
we have seen this change to affect the numerical value of the crossover temper- 
ature To in only about the second significant digit. In.our future work we want 
to focus on a numerical analysis of the HamiltonianLL3 and the QTM approach 
(to be published) yielding identical numerical values for Tq. 

Second, even within the level crossing scenario of our analytic reasoning the 
result (p8|) represents an analytic approximation to the characteristic temper- 
ature To. A comprehensive numerical treatment of ( |l5| . [l6| ) beyond the (10A) 
approximation is on the way. 

As a result of the analytical and numerical computations of finite N— 
QTM's we expect a decisive answer to the open problems. 
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